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ABSTRACT 

Long Baseline Array imaging of the z=0.663 broad line radio galaxy PKS 1421-490 reveals a 400 
pc diameter high surface brightness hotspot at a projected distance of approximately 40 kpc from the 
active galactic nucleus. The isotropic X-ray luminosity of the hotspot, £2-10 koV = 3 x 10"*^ ergs s~^, 
is comparable to the isotropic X-ray luminosity of the entire X-ray jet of PKS 0637-752, and the 
peak radio surface brightness is hundreds of times greater than that of the brightest hotspot in 
Cygnus A. We model the radio to X-ray spectral energy distribution using a one-zone synchrotron 
self Compton model with a near equipartition magnetic field strength of 3 mG. There is a strong 
brightness asymmetry between the approaching and receding hotspots and the hot spot spectrum 
remains flat (a ~ 0.5) well beyond the predicted cooling break for a 3 mC magnetic field, indicating 
that the hotspot emission may be Doppler beamed. A high plasma velocity beyond the terminal jet 
shock could be the result of a dynamically important magnetic field in the jet. There is a change 
in the slope of the hotspot radio spectrum at GHz frequencies from a ^ 0.5 to a < 0.2, which we 
model by incorporating a cut-off in the electron energy distribution at 7min « 650, with higher values 
implied if the hotspot emission is Doppler beamed. We show that a sharp decrease in the electron 
number density below a Lorentz factor of 650 would arise from the dissipation of bulk kinetic energy 
in an electron/proton jet with a Lorentz factor Fjet ^ 5. 

Subject headings: galaxies:active - galaxies:jets - quasars individual (PKS 1421-490) 



1. INTRODUCTION 

PKS 1421-490 was f irst reported as a bright, flat spec- 
trum radio source bv lEkeii (fT969h . Subsequent VLBI 
imaging revealed lOmas scale structure within the bright- 
est component of this source (jPreston et al. "l989) . Stud- 
ies at the Australia Telescope Compact Array (ATCA) 
later revealed significant radio emission on arcsecond 
scales extending s outh-west from the brightest compo- 
; nent (|Lovelllll997l ). For this reason, PKS 1421-490 was 
included in a Chandra survey of flat spectrum radio 
' quasa r s with arcsecond sca le radio jets (Mar shall et al.l 
; l2005h . iGelbord et all (|20f)5h (from here on G05) reported 
on recent X-ray ( Chandra) , optical (Magellan) and radio 
(ATCA) imaging of this source. We refer the reader to 
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that paper for the details of these observations and im- 
ages. 

Figure [T] illustrates the arcsecond scale radio structure 
of PKS 1421-490; it is annotated to show the naming 
convention for different components in the radio image 
used by G05, as well as the correct interpretation of each 
of these components brought out in this study. G05 ob- 
tained an optical spectrum of region B, and suggested it 
was not associated with an active galactic nucleus (AGN) 
in view of the apparent lack of spectral lines (due to 
poor signal to noise ratio in that spectrum). Region 
A wa s known to co ntain bright VLBI scale radio struc- 
ture (jPreston et al.l il989) and had a flat radio spectrum 
{a < 0.5). Region B was also known to be much weaker 
than region A at radio wavelengths. Consequently region 
A was thought to be an AGN, while region B was (erro- 
neously) interpreted by G05 as a jet knot. In this paper 
we show that in fact region B is the active galactic nu- 
cleus (see ®, and that region A contains a high surface 
brightness hotspot. The main focus of this paper is the 
interpretation and modeling of the exceptional hotspot 
in region A which has until now been interpreted as an 
AGN. 

One of the major results of this paper relates to an 
observed low frequency flattening in the hotspot radio 
spectrum at GHz frequencies, which indicates that the 
underlying electron energy distribution flattens towards 
lower energies. The low energy electron distribution is 
not only important for calculating parameters such as 
the number density and energy density, it also provides 
important constraints on the particle acceleration mech- 
anism. In turn, this will help to address more fundamen- 
tal issues such as jet composition and speed. We now 
present a brief overview of the literature relating to the 
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Fig. 1.— ATCA 20.2 GHz image of PKS 1421-490 with source 
components labeled. This image was first presented in G05. To 
avoid confusion, the naming convention used by G05 is also in- 
cluded. Contour levels: 1.5 mjy/beam X (1, 2, 4, 8, 16, 32, 
64, 128, 256, 512, 1024). Peak surface brightness: 2.11 Jy/beam. 
Beam FWHM: 0.54 X 0.36 arcsec. The scale of the image is 7.0 
kpc/arcsecond. 

low energy electron distribution in jets and hotspots. 

In a small number of objects, flattening of the hotspot 
radio spectra towards lower frequencies has been ob- 
served. In most cases, synchrotron self absorption 
and free-free absorption can be ruled out, and the 
observed flattening is interpreted in term s of a turn- 
over in the el ectron energy distribution fLeah v et al.l 
119891 : iCariliret al. 1988, 1991; Lazio et al. 2006). When 
modeling hotspot spectra, the turn-over in the elec- 
tron energy distribution is usually approximated by 
setting the electron number density equal to zero be- 
low a cut-off Lorentz factor 7min- In each case where 
flattening of the hotspot radio spectrum has been di- 
rectly observed, estimates of 7min sue of the order 
of several hundred: Cygnus A, 7min ~ 300 - 40 
(iCariin et al.lfT99ll:lLazio etanr2 006: Har dcastlell2001bn: 
3C295 7inin ~ 800 ([Harris e t al. 2000'; 'Hardcastl(^ 
200TD): 8C123, 7 rr,in ~ 1000 fHar^ 
Hardcastld [2001H) : PKS 1421-490 



dcastle ct al. 2001a; 
7min ^ 650 (this 

w ork) . 

iLeahv et al.l (|1989[ ) presented evidence for a low en- 
ergy cut-off in two other hotspots; 3C268.1 and 3C68.1. 
In both these ojects, the hotspot radio spectra are sig- 
nificantly flatter between ISGMHz and l.SGHz than they 
are above 1.5GHz, which suggests a similar value of 7min 
to those listed above, provided the ho tspot magn e tic field 
strengths are similar. More rece ntlv lHardcastli (|2001bD 
reported on a possible detection of an optical inverse 
Compton hotspot in the quasar 3C196. By modeling 
the synchrotron self Compton emission and assuming a 
magnetic field strength close to the equipartition value. 



they inferred a cut-off Lorentz factor 7i„in ~ 500. All of 
the above listed 7min estimates appear to be distributed 
aroun d a v alue of 7min 600. However, iBlundell et al.l 
(|2006D and lErlund et all (|2008[ ) have inferred the exis- 
tence of a low energy cut-off at 7min ^ 10^ in the hotspots 
of the giant radio galaxy 6C 0905-3955. Their method of 
detecting the low energy cut-off is quite different to those 
described above, and is based on the interpretation of an 
absence of X-ray emission from the eastern hotspot and 
radio lobe in that source. 

In f|8] we show that a turn-over in the electron energy 
distribution at 7,nin ~ 100 — 1000 can arise naturally 
from dissipation of jet energy if the jet has a high proton 
fraction and a bulk Lorentz factor Fjet ^ 5. However, 
IStawarz et al.l (|2007D have suggested that the low fre- 
quency flattening in the radio spectrum of a Cygnus A 
hotspot is not related to the turn-over in electron energy 
distribution. Rather, they argue, it indicates a transition 
between two different acceleration mechanisms. 

Electron energy distributions with a low energy cut-off 
have also been discussed in relation to pc-scale jets. The 
absence of significant Faraday depolarization in compact 
sources suggests that the number density of electrons 
with Lorentz facto r 7 > 100 grea tl y exceeds that of lowe r 
energy particles dW ardle' '1 9771: Ijones fc Odell 119771) . 
iGopal-Krishna et al . (2004) have argued that some sta- 
tistical trends in superluminal pc-scale jets may be un- 
derstood in terms of effects arising from a low energy 
cut-of f in the electron energy distribution. iTsang &: Kirtd 
(|2007( ) have suggested that a low energy cut-off in the 
electron spectrum can alleviate several theoretical diffi- 
culties associated with the inverse Compton catastrophe 
in compact radio sources, including anomalously high 
brightness temperatures and the apparent lack of clus- 
tering of powerful sources at 10^^ K. However, circular 
polarization in the pc-scale jet of 3C279 requires a min- 
imum Lorentz factor 7mi„ < 20 (jWardle et al.lll998f ). 

Observational constraints on the low energy electron 
distribution in extragalactic jets on kpc-scales are rare. 
A low energy cut-off in the electron energy distribution at 
7min ^ 20 has been estimated for the jet of PKS 0637-752 
(500 kpc from the nucleus) through modeling of the radio 
to X-ray spectral energy distribution in terms of inverse 
Compton scattering of the cosmic microwave background 
(jTavecchio et al.l 120001: [Uchivama et al.ll2005[ ). 

This paper is structured as follows: In ^we discuss 
our observations and data reduction. In ^}3|we discuss the 
active galactic nucleus - in particular the optical spec- 
trum and broad band spectral energy distribution. In f|4] 
we present the VLBI image of the northern hotspot and 
derive plasma parameters by modeling the broad band 
spectral energy distribution. In Sj5| we independently es- 
timate the hotspot plasma parameters by modeling the 
radio spectrum of the entire radio galaxy. In Sj6] we dis- 
cuss the incompatibility of the observed spectrum with 
the standard continuous injection plus synchrotron cool- 
ing model for hotspots. In SjTjwe consider Doppler beam- 
ing as a possible cause of the high radio surface bright- 
ness and various other properties of the hotspot. In fJ5] 
we consider the dissipation of energy associated with a 
cold proton/electron jet and present an expression that 
relates the energy of the peak in the electron energy dis- 
tribution to the jet bulk Lorentz factor. We then con- 
sider the implications of this expression in the case of the 
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northern hotspot of PKS 1421-490 and other objects. In 
f|9]we summarize our findings. 

Throughout this paper we assume cosmology fl^ = 
0.73,rJn.f = 0.27, i/n = 71 kms'^ Mpc~\ and we define 

so that the flux 



log(Fi/F2) 

log(l/l/i^2) 



the spectral index as a = 

density iv oc v^". 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. Summary 

We observed PKS1421-490 with the Long Basehne Ar- 
ray (LBA) at 2.3 and 8.4 GHz and with the ATCA at 2.3, 
4.8, 8.4 and 93.5 GHz. We have also made use of ATCA 
radio da ta (4.8, 8.6, 17 . 7 and 20.2 GHz) previously pub- 
lished in lGelbord et al.l (|2005f ) as well as archival 1.4 GHz 
ATCA data. We combined these data with previously 
published infra-red, optical and X-ray flux densities to 
construct radio to X-ray spectra for the northern hotspot 
and the core as well as an accurate radio spectrum of the 
entire radio galaxy. Table[T]lists the observation informa- 
tion and references for all data used in this study. Figure 
[5] presents the spectra and indicates the source of each 
data point. 

In addition to these data, we obtained an optical spec- 
trum of region B in order to confirm the classification of 
that region as an active galactic nucleus. The spectro- 
scopic observations are described in ! j2.5l 

As well as describing the observations and data reduc- 
tion steps, this section includes a description of some 
non-standard procedures that were required to construct 
the hotspot radio spectrum. Specifically, non-standard 
procedures were required to determine the hotspot flux 
density from the 8.4GHz LBA data-set due to limited 
(u, v) coverage. These non-standard procedures are de- 
scribed in ji2.2.3l Non-standard procedures were also 
used to obtain a lower limit to the hotspot flux density 
at 93.5GHz. This procedure is described in i j2.4l 

2.2. VLBI 
2.2.1. LBA Observations at 2.3 GHz 

PKS 1421-490 was observed with 6 elements (ATCA, 
Mopra, Parkes, Tidbinbilla 70m, Hobart and Ceduna) of 
the Australian Long Baseline Array (LBA) on March 23 
2006. A full 12 hour synthesis was obtained, recording 
a single 16 MHz bandwidth in both left and right hand 
circular polarization. Regular scans on a nearby phase 
calibrator, PKS 1424-418, were scheduled throughout 
the observation, as well as scans on a point like source, 
PKS 1519-273 (iLinfield et al.l[l989l) . used for gain cal- 
ibration. Unfortunately, due to hardware issues, we 
were only able to process right hand circular polariza- 
tion. However, this will not affect our results as we do 
not expect the hotspot emission to be significantly cir- 
cularly polarized. Circular p olarization in AGN seldom 
exceeds a few tenths of 1% (jRavner et al.l[2000l ). Data 
were recorded to VHS tapes using the S2 system and 
correlated using the LBA hardware correlator with 32 
channels and 2 second integration time. The data were 
correlated twice: once with the phase tracking centre lo- 
cated at the position of the radio peak in region A, and 
once with the phase tracking centre located 5 arcsec- 
onds away, at the position of the core (region B). 

The initial calibration of the visibility amplitudes was 
performed in AIPS, using the measured system temper- 
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Fig. 2. — Spectral energy distribution of the core and North- 
ern hotspot of PKS1421-490, as well as the radio spectrum of 
the entire radio galaxy. Symbols indicate the source of the data 
as follows: Filled squares, ATCA; open squares, LBA; open cir- 
cles, Magellan (taken from G05); filled triang les, 2MASS (taken 
from G05); open triangles, MOST (taken from 'Large et al.! fl981|) 
and Campbell-Wilson & Hunstead (1994)), open diamonds, Parkes 
(taken from Wills, 11975 )); filled circles, Chandra (taken from G05). 
The solid lines through the filled circles indicate the Itr range of 
X-ray spectral index permitted by the Chand ra d ata. Tips of ar- 
rows mark upper and lower limits (see section [2.41 for discussion of 
methods used to obtain ATCA limits). Symbol sizes are greater 
than or approximately equal to error bars except for those points 
where error bars have been plotted. 

atures and antenna gains. We obtained simultaneous 
ATCA data during our observation, and this allowed us 
to bootstrap the LBA flux scale to the ATCA flux scale 
by comparing simultaneous measurements of the point 
like source PKS 1519-273. After scaling the gains us- 
ing this bootstrapping method, and correcting the resid- 
ual delays and rates via fringe fitting, the data-set from 
the phase refer ence source was exported to DIFMAP 
(jShepherdji 19971 ) where it was edited and imaged. Am- 
plitude and phase self calibration corrections obtained 
from imaging the phase reference source were imported 
into AIPS using the cordump^^ patch kindly supplied to 
us by Emil Lenc. These phase and amplitude correc- 
tions were then applied to PKS 1421-490, and the data 
exported to DIFMAP for deconvolution and self calibra- 
tion. The resulting image is shown in figure [3l We mea- 
sure a hotspot flux density of 4.25 ± 0.2 Jy at 2.3 GHz. 
Preston et al. (1989) obtained a flux density of 4.1 Jy 
at 2.3 GHz for the northern hotspot by model fitting 
SHEVE (Southern Hemisphere VLBI Experiment) data 

12 rpj^g cordump patch i s available for DIFMAP at 
|http : //astro nomy . swin. edu. au/$\sim$ele nc/DlfmapPatches/| 
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Note. — The uncertainties in ATCA flux density are dominated by the uncertainty in the absolute flux calibration, which is estimated to be 2%, except 
for 93.5 GHz, where the uncertainty is estimated to be 10%. The lower limit on the hotspot flux at 93.5 GHz com es from mak i ng an assu mption about the 
non-hotspot spectrum extrapo l ated to higher frequencie s from 8.4 GHz . See section 12.41 for details. References: (1) iLarge et al.l (flQSTh . (2) Iwillg (flOTSh . (3) 
ICampbell-Wilson fc HunsteadI (fT99l) , (4) This work, (5) IGelbord erall (|2005l) . 

^ Convert to linear resolution using 7.0 kpc/arcsecond^ Only short baselines on which the radio galaxy is unresolved were used to measure the flux density.^ Res- 
olution of shortest baseline.*^ See scction |2.4l 
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with the simplest model consistent with the data (two 
circular Gaussians). 

We attempted to detect compact structure within re- 
gion B (the core) using the data-set that had been corre- 
lated with the phase centre at that position. The time- 
averaging- and bandwidth-smeared emission from region 
A in this dataset was first cleaned to remove the side- 
lobes, but we were unable to detect any emission from 
the location of the core to a limit of approximately 8mJy 
(5 fj). The upper limit to the flux den sity of region B at 
4.8 GHz is 7 mJy (iGelbord et al.l[200l 

2.2.2. LB A Observations at 8.4 GHz 

PKS 1421-490 was observed with 5 elements of the 
LBA (ATCA, Mopra, Parkcs, Hobart and Ceduna) in 
May 2005 at 8.4 GHz. A fuU 12 hour synthesis was ob- 
tained, recording a single 16 MHz bandwidth in both left 
and right hand circular polarization. At this frequency, 
the Northern hotspot of PKS 1421-490 is completely re- 
solved on all but the shortest three baselines (baselines 
between Parkes, the ATCA and Mopra). We fringe fit- 
ted our target source data using a point source model 
in AIPS, before performing model-fitting and phase self 
calibration iterations in DIFMAP. 

2.2.3. Determination of Hotspot Flux Density from the 
8.4 GHz LBA Data Set 

Due to the small number of baselines, we use model- 
fitting in the {u, ?;)-plane rather than CLEAN deconvo- 
lution to measure the hotspot flux density at 8.4 GHz. 
Model-fitting involves specifying a starting model in the 
image plane, consisting of a number of elliptical Gaussian 
components, each with a particular flux density, position, 
size and position angle, then allowing the model fitting 
algorithm to locate a chi-squared minimum by fitting the 
Fourier transform of the model to the {u, 'i;)-data. 

Care is necessary when comparing the LBA flux den- 
sity measurements at the two different frequencies, due 
to the limited {u, v) coverage. At 2.3 GHz, the data 
cover {u, v) spacings between 0.5 and 13 MA, while at 
8.4 GHz the data cover {u, v) spacings between 2 and 9 
MA. Therefore, provided the source structure can be de- 
scribed by a simple model consisting of a set of Gaussian 
components, the comparison of model flux densities will 
be vahd. 

In order to determine the range of allowable flux densi- 
ties in the 8.4 GHz data-set, we specified a wide range of 
different models consisting of 3, 4 or 5 elliptical Gaussian 
components, broadly consistent with the 2.3 GHz image, 
then let the model-fitting algorithm adjust the model to 
fit the 8.4 GHz data. While it is not possible to pre- 
cisely constrain the flux density of the hotspot with only 
three baselines, we found that the total flux density of 
all acceptable models (using between three and five el- 
liptical Gaussian components and a wide range of initial 
model parameters) was never less than 2.9 Jy, and the 
flux density of the best-fitting model was 3.2 Jy. The 
8.6 GHz ATCA image contains an unresolved source of 
3.3 Jy at the position of the hotspot, and this provides 
an upper limit to the hotspot flux density at 8.6 GHz. 
We therefore ado|)t a hotspot flux density at 8.4 GHz of 
^8.4GHz = 3.2to;3 Jy. 



2.3. ATCA Observations 

ATCA observations of PKS 1421-490 were made simul- 
taneously during our LBA observations. We recorded 
a single 64 MHz bandwidth at 8.4 GHz and 128 MHz 
bandwidth at 4.8 GHz during our first LBA observation 
in 2005. A single 128 MHz bandwidth at 2.3 GHz was 
recorded during our second LBA observation in 2006. For 
each of these observations, the ATCA was in a compact 
configuration so we could not image the source in detail, 
but we were able to obtain accurate total source flux 
density measurements (see Table [T]). Total source flux 
density measurements were also obtained at 1.4 GHz us- 
ing archival ATCA data. Standard calibration and imag- 
ing procedures were used with the MIRIAD processing 
software. 

In August 2005 we obtained a full 12 hour synthe- 
sis with the new 3mm receivers. Again, standard cali- 
bration and imaging procedures were used in MIRIAD. 
The flux density scale was determined from scans on the 
planet Uranus and confirmed using the point-like source 
PKS 1921-293, the flux density of which had been mea- 
sured 4 days prior to our observing run as 8.8 ± 0.9 Jy 
at 93.5 GHz. We detected a single point like component 
in the 93.5GHz image of PKS 1421-490 coincident with 
region A, the flux density of which we regard as being 
the total source flux density at this frequency, since the 
resolution of the shortest baseline larger than the source. 
The upper limit on flux density at 93.5 GHz for region B 
and C is 5 mJy (5 a). 

Errors in the flux densities reported in Table [1] are 
dominated by uncertainties in the primary flux calibra- 
tion which are estimated to be of order ~2% at cm wave- 
lengths, and of order ^10% at 3mm. 

2.4. Constraints on the Hotspot Radio Spectrum from 
ATCA Images 

The hotspot is unresolved in the ATCA images, and is 
blended with emission from the surrounding regions. We 
are therefore unable to directly measure the flux density 
of the hotspot from the ATCA data. However, we are 
able to constrain the flux density of the hotspot, and we 
now discuss the methods used to obtain upper and lower 
limits. 

Using radio dat a that was first presented in 
IGelbord et all ((20051) . we find upper limits on the hotspot 
fiux density at 17.7 GHz and 20.2 GHz by summing the 
CLEAN components at the position of the hotspot. Simi- 
larly, we obtain an upper limit to the hotspot flux density 
at 93.5 GHz from the measurement of total source flux 
density at that frequency. We obtain a lower limit on the 
hotspot flux density at 93.5 GHz via the following steps: 

1. We subtract the LB A- measured hotspot flux den- 
sity from the total source flux density at 2.3 GHz 
and 8.4 GHz to obtain two estimates of the non- 
hotspot flux density, from which we calculate a 
non-hotspot spectral index (aj'g Qy^(^o^"^^otspot) 
= 0.78). 

2. We extrapolate the non-hotspot power law to 
93.5 GHz. 

3. We reasonably assume that the non-hotspot spec- 
trum becomes steeper towards higher frequencies. 
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Fig. 3.— VLSI image of the Northern Hotspot of PKS 1421- 
490 at 2.3 GHz. Contour levels: 3.5 mjy/beam X (1, 2, 4, 8, 
16, 32, 64, 128). Peak surface brightness: 0.47 Jy/beam. Beam 
FWHM: 13.5mas X 11.6mas. The scale of this image is 7.0 pc/milli- 
arcsecond. 



Therefore the extrapolated flux density from step 
2 is an upper limit to the non-hotspot flux density. 

4. We subtract the non-hotspot upper limit from the 
observed entire source flux density to obtain a lower 
limit on the hotspot flux density at 93.5 GHz. 

The assumption in step 3 is based on the observation that 
the non-hotspot emission arises in the lobes, jets and the 
southern hotspot (the core is neghgible). Jet and lobe 
spectra are often observed to steepen towards higher fre- 
quency. Indeed, the 17.7 GHz and 20.2 GHz ATCA im- 
ages indicate that the spectral index of the northern lobe 
region steepens signiflcantly at higher frequency. The 
limits on hotspot flux density obtained from the ATCA 
images are represented by the tips of the arrows in Figure 
[21 and are listed in Table [TJ 

2.5. Optical Spectroscopy 

Optical spectra were taken with the Magellan IMACS 
camera on 14 May 2005 in service mode. Three ten 
minute exposures were obtained using a long slit (0'.'9 
width) aligned with regions A and B. The 300 lines/mm 
grism was used, to yield a spectral resolution of i? ^ 1000 
spanning roughly 4000-10000 A. The spectra were re- 
duced with IRAF. No standard stars were observed, so 
no effort was made to flux calibrate the spectra or to 
remove telluric absorption features. 

No signiflcant spectral features were detected in the 
spectrum of region A, consistent with synchrotron emis- 
sion from a hotspot. However, only a very high equiva- 
lent width emission line could have been detected due to 
the low signal to noise ratio of these data. 

The spectrum of region B contains several broad and 
narrow emission lines, which allowed a precise determi- 
nation of the redshift (see §3.1|1 . This is not the first 
spectrum of the nucleus to be published — a spectrum 
of region B was presented in G05. However, the spec- 
trum presented in G05 did not allow identification of any 
spectral lines because the strongest spectral features fell 
beyond the wavelength coverage, and the spectrum was 
taken as a single short exposure with a high background 
due to the pre-dawn sky. 
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Fig. 4. — The normalized optical spectrum of region B. This spec- 
trum clearly shows that region B is a broad line AGN at redshift 
0.6628. Increasing in wavelength, the arrows indicate the positions 
of the following emission lines: Mg II 2799, [Ne V] 3426, [O II] 
3726/3729, [Ne III] 3869, H/3, [O III] 4959, [O III] 5007. 

3. REGION B: THE ACTIVE GALACTIC NUCLEUS 

3.1. Optical Spectrum oj Region B 

The normalized optical spectrum of region B is dis- 
played in Figure [H We detect several broad and nar- 
row emission lines: Mg II 2799, [Ne V] 3346 and 3426, 
the blended [O II] 3726,3729 doublet, [Ne III] 3869 and 
3967, H(5 (marginal), H7, [O HI] 4363, H/3, and [O HI] 
4959 and 5007. H/3 has both a broad and narrow com- 
ponent; their measured FWHM values (uncorrected for 
instrumental resolution) are 6500 km/s and 516 km/s re- 
spectively. The narrow H/3 line width is consistent with 
that of the OIII lines (510 km/s). From these features 
we measure the redshift z — 0.6628 ± 0.0001. 

3.2. Spectral Energy Distribution of the Core 

We did not detect the active galactic nucleus with the 
LBA, but obtain an upper limit on the fiux density of 
approximately BmJy at 2.3 GHz. This is consistent with 
the ATCA core fiux density measurements (see Table [Ij . 
The core is completely dominated by its optical emis- 
sion. In fact, the optical fiux density is so great relative 
to the radio {afl^^^ = 0.23 ±0.02), the core would 
be classified as radio quiet in the strictest sense. The 
radio to optical spectral index cannot be explained in 
terms of standard sync hrotron self abs orption models for 
fiat radio spectra fee. iMarscheij|1988l ). since this would 
require at least part of the jet to be self absorbed at opti- 
cal wavelengths and would imply an unrealistically high 
magnetic field strength. This region has an optical spec- 
tral index ctp = 0.2 z fc 0.1, in the range typical of quasars 
(|Francis et al.l 119911 ). suggesting that the strong optical 
emission may be due to an unusually large contribution 
from the accretion disk thermal component. The radio 
to X-ray spectral index (q;^|'=^hz = O-^IO ± 0.005) is typ- 



ical o f radio galaxies at similar redshift feg. lBelsole et al.l 
l2006f ). The optical to X-ray spectral index is Uox = 1-62 
(G05). There is clearly an excess of optical fiux relative 
to the radio and X-ray fl ux when compared t o samples of 
other radio galaxies feg. lGambill et ani2003f) . Note that 
measurements of the B-band magnitude have shown no 
variability, to within 0.6 magnitudes, over the past 35 
years (G05). 

4. THE NORTHERN HOTSPOT 

4.1. Morphology 

Figure [3] shows the LBA image of the Northern hotspot 
at 2.3 GHz. Less than 0.2 Jy (5% of the hotspot flux den- 
sity) remains on the longest baseline (~ 12.9MA), imply- 
ing that there is little structure on scales smaller than 
15 mas (100 pc). This limit on substructure within the 
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hotspot is relevant to possible synchrotron self absorp- 
tion models for the hotspot spectrum, which we discuss 
further in gUl 

The flux density of the hotspot in our LBA image 
is 60 % of the total source flux density at 2.3 GHz, 
and 75 % at 8 GHz. The peak surface brightness is 
-^2 Tghz = 2600 Jy/arcsec^. Extrapolating to 8 GHz as- 
suming Ii, (X iy~'^'^ (the spectral index of a ~ 0.2 between 



these frequencies is calculated in H4.2p and accounting 
for cosmological dimming and redshifting, we find that 
the peak surface brightness of the northern hotspot of 
PKS 1421-490 would be more than 1000 times brighter 
than the brightest hots pot of Cygnus A if they were at 
the same redshift (C arilli et al.l 119991 ). The monochro- 
matic hotspot luminosity is L2.3GHZ = 8 x 10^'' WHz~^. 

A protrusion on the Eastern edge of the hotspot re- 
sembles the "co mpact protrus ions" seen in n umerical 
simula tions (eg. iNormanI |1996[ ) . According to iNormanI 
()1996f ). a compact protrusion is produced in their 3- 
D non-relativistic hydrodynamic simulations when the 
light, supersonic jet reaches the leading contact discon- 
tinuity. At this point, the jet is generally flattened to a 
width substantially less than the inlet jet diameter, and 
the compact protrusion arises where the jet impinges on 
the contact discontinuity surface. 

The width of the hotspot at the peak (region HSl) is 
400pc measured perpendicular to the inferred jet direc- 
tion. The lower surface brightness emission behind the 
hotspot peak (region HS2) is 700pc measured perpen- 
dicular to the jet direction. The length of the hotspot 
(regions HSl and HS2) is approximately Ikpc. The geo- 
metri c mean of the maj o r and minor axes (for comparison 
with Har dcastle et al.l (|1998 fl and 'Jevakum ar fc Saikial 
diOOO)) is0.63kpc. This is a factor of 4 below the median 
value (2.4 kpc) of hotspot sizes given in Hardcastle et al.l 
(|1998f ). However, the size of the hotspot relative to the 
linear size of the source is consistent with the correla- 
tion betwee n these p arameters given in lHardcastle et all 
([l99S) and Jeyakumar & Saikia ( 206^. 

The jet exhibits a bend of almost 60 degrees (pro- 
jected) approximately 5 arcseconds (35kpc) from the core 
at the western end of the ridg e of emission extend ing west 
from the hotspot in Figure [TJ iBridle et all (|1994D showed 
that hotspot brightness is anti-correlated with apparent 
jet deflection angle. They found that, for the twelve 
quasars in their sample, the ratio of hotspot flux den- 
sity to lobe flux density decreases with larger jet bending 
angles, particularly when the deflection occurs abruptly. 
PKS 1421-490 does not follow this trend. 

We detect what appears to be a jet knot (region Jl) 
512 mas (~ 3.5 kpc projected) at position angle -107 de- 
grees (North through East) from the hotspot peak. The 
knot is extended along a position angle almost perpen- 
dicular to the apparent jet direction. The major axis of 
the knot is poorly constrained due to the low signal to 
noise of this component, but the data suggests a width 
of approximately 400 - 600pc. 

4.1.1. Interpretation of Region HS2 

We now consider the interpretation of the lower surface 
brightness region HS2 just behind the hotspot peak. As 
mentioned above, the diameter of region HS2 perpendic- 
ular to the jet direction (700pc) is much larger than the 
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Fig. 5. — Radio spectrum of the entire radio galaxy (filled 
squares) and hotspot (filled circles and tips of arrows). This figure 
serves to illustrate the low frequency flattening in the hotspot radio 
spectrum referred to in the text f i|4.2.1l and ^El . The dotted lines 
illustrate power law fits to specific sections of the data. 

diameter of region HSl (400pc). The surface brightness 
of region HS2 is more than a factor of 10 times the peak 
surface brightness of the brightest hotspot of Cygnus A. 
In addition, the flux density from region HS2 alone 
0.6 Jy at 2.3 GHz) is more than 4 times the total flux 
density of the whole counter lobe and hotspot. There 
are two possible interpretations for region HS2, and the 
interpretation of this region has implications for the in- 
terpretation of region HSl. 

The first interpretation is in terms of emission from 
turbulent back-flow in the cocoon. If this interpretation 
is correct, we cannot appeal to Doppler beaming to ac- 
count for the high surface brightness of region HS2 rela- 
tive to other hotspots, and the high flux density relative 
to the counter hotspot and lobe. If we cannot appeal 
to Doppler beaming for region HS2, it would seem un- 
reasonable to appeal to Doppler beaming to explain the 
high surface brightness of region HSl. The arm length 
symmetry places a tight upper limit on the expansion 
velocity of the lobes at ^expansion < 0.1c, indicating that 
the whole complex (region HSl and HS2) cannot be ad- 
vancing relativistically. 

The second possible interpretation for region HS2 is 
that the emission is associated with oblique shocks in the 
jet as it approaches the hotspot. In this case, we may 
appeal to Doppler beaming to explain the high surface 
brightness of both regions HSl and HS2. However, this 
interpretation would imply that the jet diameter at HS2 
(~700pc) is signiflcantly greater than the diameter of the 
hotspot at HSl (~400pc) and also greater than the jet 
diameter at Jl (~ 400 — 600pc). It should be noted that 
the width of region Jl, presumably associated with a jet 
knot, is poorly constrained due to the low signal to noise 
of this component. Future LBA observations at 1.4GHz 
may provide better constraints on the size of regions Jl 
and HS2. 

4.2. Modeling the Hotspot Spectral Energy Distribution 
4.2.1. Low Frequency Flattening 

Figure [5] illustrates that the hotspot radio spectrum 
changes slope at GHz frequencies, becoming flatter to- 
wards lower frequency. We now discuss this feature in 
more detail and consider the possible causes. 

The hotspot spectral index calculated from our LBA 
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flux density measurements is relatively flat at GHz fre- 
quencies (ai:|§Hz = 0.22l[J[J^). The hotspot spec- 
trum cannot continue with this slope to millimeter wave- 
lengths, since it would substantially over-predict the ob- 
served 93.5 GHz flux density. We therefore require that 
the hotspot spectrum be steeper at frequencies above 
8 GHz , with spectral index a > atiiA^AGiii^'" = 0-48- 

Our conclusion of a flat spectral index at GHz fre- 
quencies based on the LBA flux density measurements is 
strengthened by inspection of the whole source spectrum 
(see Figure[5]). The flux density of the entire radio galaxy 
has a spectral index of 0.58 above 8.4 GHz, but flattens 
to a spectral index of 0.34 below 4.8 GHz. The northern 
hotspot is the dominant component at GHz frequencies, 
hence, the flattening of the total source spectrum implies 
there is flattening in the hotspot spectrum. There are a 
number of possible causes of this GHz frequency flatten- 
ing, but most of them are implausible. We now consider 
a number of such explanations. 

If synchrotron self absorption were responsible for the 
flattening, the required magnetic field strength is Bq ^ 
lO~^v^e^Fp^{l + where Bq is the magnetic field 
strength in Gauss, Fp is the peak flux density in Jy, 
is the frequency in GH z at the peak an d is the angular 
size in milliarcseconds (de Yoimg 200^, pg. 325). In the 
case of the Northern hotspot of PKS 1421-490 we esti- 
mate (conservatively) Vp ^ IGHz, 6 ~ lOOmas, Fp ~ 5Jy 
and z = 0.663. Therefore, a magnetic field strength of 
i? ^ 20 G is required to produce the observed flattening 
— four orders of magnitude greater than the equipar- 
tition magnetic field strength. Less than 0.2Jy (5% of 
the hotspot flux density) remains on the longest base- 
lines (12.9 MA ^15 mas resolution), implying that the 
hotspot cannot be composed of many small self-absorbed 
sub-components. Therefore, we do not consider syn- 
chrotron self absorption to be a viable explanation for 
the flattening. 

We next consider free-free absorption by interstellar 
clouds in the hotspot environment as a possible mech- 
anism for the observed flattening of the radio spec- 
trum. Consider a cloud of size ikpc kpc, temperature 
T4 X 10'^ K, electron number density Ue cm~^ and pres- 
sure p-12 X 10~^^ dyn/cm^. The optical depth t to free- 
free absorption at a frequency i^ghz GHz is given by (eg. 
Ide Yound^oM pg.326) 

r = 3.3 X lO-^n^LkpcJ^GHz T^4"'-'' (1) 

= 2xlO-^Lkpc.GHzV^2 i^yT,^'-'' (2) 

Assuming a characteristic pressure in the outer regions 
of an elliptical galaxy p ~ 10"^^ dyn/cm^, a character- 
istic temperature for an ionized cloud T ~ 10** K and 
a reasonable cloud size L < 1 kpc, the optical depth to 
free- free absorption above 1 GHz is less than 5 x 10~^. 

We therefore interpret the change of slope in the 
hotspot radio spectrum in terms of a change in the un- 
derlying electron energy distribution. In i i4.2.4l wc model 
the SED by incorporating a low energy cut-off in the elec- 
tron energy distribution. A low-energy cut-off at some 
minimum Lorentz factor 7min produces a spectrum with 
Fi, (X v^^^ at frequencies below the characteristic emis- 
sion fregu£n£;Y_o£elBctrons with Lorentz factor 7min (see 
eg. IWorrall fc Birkinshawl l2006f ) . We must emphasize 



that an instantaneous cut-off in number density is not 
physical - it is merely an approximation to a sharp turn- 
over in the electron energy distribution. In SJ5] we show 
that interpreting the observed flattening in terms of a 
turn-over in the electron energy distribution has consid- 
erable implications. 

Low frequency flattening in hotspot spectra has been 
observed in a small number of other objects (see flj- 

4.2.2. The High Frequency Synchrotron Spectrum 

The hotspot spectrum remains relatively flat between 
8GHz and 93.5GHz, having spectral index 0.4 < a < 0.6 
(based on the two point spectral index from the 8.4GHz 
LBA data point to the ATCA upper and lower limits at 
93.5GHz). The synchrotron spectrum above 93 GHz is 
poorly constrained, but the simplest model — a power 
law spectrum with spectral index 0.4 < a < 0.6 and an 
exponential cut-off at high frequency (i.e. a synchrotron 
spectrum from a power law electron energy distribution 
with number density set to zero above 7rnax) , is unable to 
satisfy the optical data and the 2MASS infra-red upper 
limits simultaneously. Either a break to a steeper spec- 
trum somewhere between ~ 10^^ — 10^'^ Hz is required, 
or a gradual cut-off at high electron energy, rather than 
an abrupt cut-off at 7max, must exist. Given the high ra- 
dio luminosity, hence high magnetic field strength, syn- 
chrotron losses are likely to be important. We there- 
fore allow for a synchrotron cooling break at an arbitrary 
break frequency in order to fit the radio through optical 
spectrum. Different choices of model spectrum are pos- 
sible, but they would not significantly affect our major 
results. We discuss the model electron energy distribu- 
tion in i )4.2.41 and further discuss the self-consistency of 
this model in fj6l 

4.2.3. Hotspot X-ray Emission 

Figure [6] shows the spectral energy distribution (SED) 
of the northern hotspot. The level of X-ray flux density 
relative to the optical flux density indicates the pres- 
ence of two distinct spectral components: synchrotron 
emission from radio to optical frequencies, and inverse 
Compton emission at X-ray frequencies and above. 

The energy density of the locally generated syn- 
chrotron emission within the hotspot (assuming no 
Doppler beaming) is more than 10'* times the energy den- 
sity of the cosmic microwave background (CMB) at this 
redshift. If the hotspot plasma is moving relativistically 
with velocity u = /3c at an angle 9 to the line of sight, the 
ratio of synchrotron to CMB energy density is reduced 
by a factor of ~ p-^^-^ ^^lere S = [r(l - Pcos0)]-^ 
is the Doppler factor and T = (f — /3^)^*/^ is the 
bulk Lorentz factor, and we have assumed the hotspot 
is associated with plasma moving through a station- 
ary volume /pattern rather than a m oving blob, so that 
= S^Il, (jLind fc Blandfordl[l985l ). Therefore, if 5 < 6 
(assuming F ~ i5), inverse Compton scattering of locally 
generated synchrotron photons is the dominant source 
of inverse Compton X-ray emission. While a Lorentz 
factor of F ^ 6 is not ruled out, such a high Lorentz 
factor is not required by the data, and we consider only 
Lorentz factors F < 3. We therefore ignore the inverse 
Compton scattering of CMB photons in the following 
treatment. We also ignore any contribution to the X- 
ray flux density from "upstream Compton" scattering 
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Fig. 6. — Hotspot Spectral Energy Distribution with SSC model 
overlaid. The "bow-tie" around the X-ray data point indicates the 
1(T range of X-ray slopes. The solid line is the best fit synchrotron 
plus self Compton model spectrum with the Doppler factor fixed 
at (S = 1. The dashed line is the best fit SSC component with 
the Doppler factor fixed at 5 = 3. The model synchrotron spectra 
for 5 = 3 and <5 = 1 are exactly the same, and so do not appear 
as separate curves in the plot. The LBA data points (plotted as 
filled circles) have error bars smaller than the symbol size. Tips of 
arrows mark the position of upper and lower limits . X-ra y, optical 
and infra-red points are taken from lGelbord et al.l | |2005|) . 

(jGeorganopoulos h Kazanail2003D . whereby electrons in 
the jet upstream from the hotspot inverse Compton scat- 
ter synchrotron photons produced within the hotspot. 
We note that if the upstream Compton process makes a 
significant contribution to the observed X-ray flux den- 
sity, the SSC flux density must be less than the observed 
flux density, in which case the magnetic field strength in 
the hotspot would be greater than that reported in Table 
[21 and therefore greater than the equipartition value. 

4.2.4. Synchrotron Self Compton Modeling 

To model the radio to X-ray spectral energy distribu- 
tion we use the standard one-zone SSC model: a spheri- 
cal region of plasma with uniform density and magnetic 
field strength. We assume that the magnetic field is "tan- 
gled" with an isotropic distribution of field direction. We 
further assume that the number density of electrons per 
unit Lorentz factor is described by 



iV(7) = 
where 




N{'~f) is the volume averaged energy distribution pro- 
duced by continuous injection of a power-law energy dis- 
tribution Af(7) — i^eT"" at a shock with synchrotron 
cooling in a uniform magnetic field downstream. It de- 
scribes a broken power-law spectrum with the electron 
spectral index smoothly changing from -a to -(a-|-l) at 
7 « 7b. The break in the electron spectrum at 7b 
corresponds to the electron energy at which the syn- 
chrotron cooling timescale is comparable to the dynam- 
ical timescale for electrons to escape from the hotspot. 
The synchrotron cooling break is discussed further in ^J6l 
For the electron energy distribution described by N{^), 
and given a particular radius, redshift, Doppler factor 



and spectral index, the synchrotron plus self Comp- 
ton spectrum is characterized by the five parameters 
Ke, B, 7min, 76, 7max- In Calculating the model spectrum, 
these parameters appear in the following combinations 
(see appendix [X]): 



A 

^syn 




(5) 


A 

^ssc 




(6) 




~ (l + 2)47r 


(7) 




= (l + z)47r""^'' 


(8) 


V2 


~ (1 + z) 47r "'^'"^'^ 


(9) 



where a = (a — l)/2 is the radio spectral index between 
frequencies vi and Vb, is the non-relativistic gyro- 
frequency. The parameters ^i, Vb and 1^2 correspond to 
the characteristic frequency emitted by electrons with 
Lorentz factor 7inin, 7b and 7niax in a magnetic field of 
fiux density B, and are therefore identified with the low 
frequency turn-over, synchrotron cooling break and high 
frequency cut-off respectively. The advantage of this for- 
mulation, described in Appendix]^ is that it allows the 
model to be specified in terms of the observed values of 
vi^Vb and U2. The parameters Agyn and Agsc are nor- 
malization factors for the synchrotron and SSC spectral 
components respectively. 

We estimate best fit values for B, K^, 7min, 7b and 7max 
using chi-squared minimization with the following three 
constraints: (1) We fix the electron energy index at a = 
2.06 so that the spectral index a = 0.53 for frequencies 
I'l << ly << i^b (that is, between about 10 GHz and 100 
GHz) . This is the electron energy index determined from 
modeling the hotspot radio spectrum as described in SJS) 
The spectral index a — 0.53 also agrees with the ratio of 
peak surface brightness (at the location of the hotspot) in 
the 17.7GHz and 20.2 GHz ATCA images. (2) We fix the 
radius at R 320pc. The radio hotspot is elongated in an 
approximately cylindrical shape of volume y « 4 x 10^^ 
m'^. A radius of 320pc gives a spherical model of equal 
volume. (3) We fix the break frequency at Vb — 500 GHz. 
This is close to the lowest break frequency allowed by the 
data. Higher break frequencies are permitted but cause 
a worse fit to the optical data. The break frequency is 
not well constrained by the data, but the results are not 
sensitive to the assumed value of the break frequency. 
(4) We fix the upper cut-off frequency at t'max = lO^** Hz 
to fit the optical fiux densities. 

We determined best fit parameter values while fixing 
the Doppler factor at i5 = 1,2 and 3. The derived model 
parameters are presented in Table [2] The uncertainties 
in Table [2] are determined from the range of parameter 
values in the set of models having < Xmin + 2.71. The 
observed X-ray spectral index was not included in the 
chi-squared calculations, but the model X-ray spectral 
index is consistent with the observed value within the 
uncertainties. 

In Figure [6] we plot the observed hotspot fiux densi- 
ties with the best fit model spectra (for Doppler factors 
fixed at S — I and 6 — 3) overlaid. The simple one-zone 
model with a near equipartition magnetic field strength 
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TABLE 2 

Synchrotron Self Compton Model Parameters for Northern Hotspot 



Fixed Parameters Derived Parameters 



S K a !/(, I^niax B B/Beq rtg 7min 7b 7max a, 

[pc] [Hz] [Hz] [mG] [xlO-''^cm-3] [xlO"'] [xlO^] 



1.0 320 0.53 5 x 10^^ 10^'' 3.2 ±0.5 l-S^oi 2.7 ±0.4 650 ± 80 > 0.8 1.1 ±0.1 0.5 ± 0.1 

2.0 320 0.53 5 x 10" lO^* 1.0 ±0.2 0.75 ±0.15 2.2 ± 0.3 800 ± 100 >1 1.4 ±0.2 0.45 ±0.1 
3.0 320 0.53 5 x 10" 10^" 0.5 ±0.1 0.5 ±0.1 1.9 ±0.3 950 ± 100 >1.1 1.6 ±0.2 0.45 ±0.1 



Note. — The quoted uncertainties on model parameters arc an estimate of the level of uncertainty from model fitting, and 
correspond to the range of parameter values in the set of models having < Xmin ± 2.71. 



provides a good description of the available data. Hard- 
castle et al. (2002) used more complicated spectral and 
spatial models for three sources, and found that this did 
not have a significant effect on the derived plasma pa- 
rameters. We are therefore confident in our parameter 
estimates using this "first-order" one-zone model. 

5. MODELING THE RADIO SPECTRUM OF THE ENTIRE 
RADIO GALAXY 

We now describe a consistency check for the model of 
the hotspot radio spectrum in terms of a cut-off in the 
electron energy distribution at 7min- This check is based 
on the observed flattening in the spectrum of the entire 
radio galaxy (Figure [7]). In order to test whether the 
observed flattening is consistent with the inferred low 
energy cut-off in the electron distribution, we fit a sim- 
ple two-component model to the radio galaxy spectrum 
between 408]VIHz and 93.5 GHz. The model components 
are: (1) The synchrotron spectrum produced by the elec- 
tron energy distribution of equation ([3]) . This component 
describes emission from the hotspot. (2) A pure power- 
law approximating emission from the rest of the source. 
Note that the core flux density is negligible compared 
with that of the jets, lobes and hotspots, so that no com- 
ponent is included to represent emission from the AGN. 

The synchrotron spectrum of component 1 is calcu- 
lated using equation (jASp . We assume the same source 
volume as in ji4.2.41 fix the lab frame break frequency at 
Vh = 500 GHz and the lab frame high frequency cut-off at 
z/2 = 1-1 X lO^^Hz, consistent with the values determined 
from modeling the hotspot spectrum in M.2.41 With 
these assumed values, the parameters Vb and 1/2 do not 
affect the shape of the spectrum below about 100 GHz, 
but they weakly affect the calculation of equipartition 
magnetic field strength. The spectrum of component 1 
is therefore determined by the spectral index ai, the 
turn-over frequency vi and the synchrotron amplitude 
^syn — Ke^o^^^- Thc flux dcnsity of the second compo- 
nent is of the form 

F.,2=F2.3GHz,2(^;^^) (10) 

Chi-squared minimization was used to determine the best 
fit values for the parameters ai, Agyn, Vmin, ^2.3GHz,2 
and a2 ■ The resulting model is shown in Figure [71 This 
simple two-component model provides an excellent fit to 
the radio galaxy spectrum. Component 1 (describing 
emission from the hotspot) is in good agreement with the 
LBA flux density measurements at 2.3 GHz and 8.4 GHz. 
We emphasize that the LBA flux density measurements 




v[Hz] 

Fig. 7. — Radio spectrum of the radio galaxy PKS 1421-490. 
Open diamonds represent flux densities of the entire radio galaxy 
from the Molonglo Synthesis Radio Telescope (MOST). Filled dia- 
monds represent flux densities of the entire radio galaxy from the 
Parkes telescope. Filled circles represent flux densities of the entire 
radio galaxy from the ATCA. Open squares represent flux densities 
of thc Northern hotspot from thc LBA (see Table |3- Also shown 
is thc spectral decomposition described in section^ Component 
1 (curved dashed line) is the angle-averaged synchrotron spectrum 
from a power-law electron distribution of the form A'^(7) oc -y~2.06 
with a low energy cut-ofi^ at a Lorentz factor corresponding to 
ui = 2.6 GHz (see equation l[7]l). This component describes emis- 
sion from the hotspot, and is consistent with the LBA flux densities 
plotted as open squares. Component 2 (straight dotted line) is a 
pure power-law with spectral index 0.64. This component is an 
approximation to the emission from the rest of the source. The 
solid line is the sum of components 1 and 2. Model parameters are 
given in Table [3] 

were not included in the fitting process, but are included 
in Figure [7] for comparison with the spectrum of compo- 
nent 1. 

Again, we point out that the cut-off in number density 
per unit Lorentz factor below 7i„in used in our modeling 
is not physically realistic, it is merely an approximation 
to an electron energy distribution with a sharp turn-over. 
The success of this simple model in simultaneously ac- 
counting for the flattening in both the hotspot spectrum 
and the radio galaxy spectrum supports an interpreta- 
tion of the flattening in terms of a sharp turn-over in the 
hotspot electron energy distribution at a Lorentz factor 
of order 7 ~ 600. 

6. SYNCHROTRON COOLING BREAK 

The magnetic field strength inferred from spectral 
modeling in section 14.2.41 implies that there should be a 
synchrotron cooling break in the hotspot radio spectrum 
at ^ 1 GHz if the electron energy distribution injected at 
the shock is a pure power law. This is inconsistent with 
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TABLE 3 

Model Parameters from Fitting the Radio Spectrum of the Entire Radio Galaxy 





Component 1 (Hotspot) 


Component 2 


F(2.3 GHz ) 

[Jy] 


[GHz] [Gauss] 


F(2.3 GHz ) a 

[Jy] 


4.15 ± 0.6 


0.53t»;i 2.%+_°o% 2.05l«;f X 10-^ bbOtfoo 


3.0 ± 0.6 0.65 ±0.1 



Note. — The quoted uncertainties on model parameters are an estimate of tire level of uncertainty 
from modclfitting, and correspond to the range of parameter values in the set of models having < 
X™,„ +2.71. 



the lower Hmit on the break frequency estimated from 
spectral modeling, v\j > 500 GHz. In this section we con- 
sider the production of the cooling break, and possible 
reasons for the inconsistency. 

The standard continuous in jection hotspot model 
(|Heavens fc IVIeisenheimeil I1987D predicts that the ra- 
dio spectrum will steepen from ain to ain + 0.5 at a 
frequency, i>b, corresponding to the electron energy at 
which the synchrotron cooling time-scale TcqoI is equal 
to the dynamical time-scale Tosc for electrons to escape 
the hotspot. The break frequency is an important con- 
straint on the physics of the hotspot. In general, it de- 
pends on the magnetic field strength, hotspot radius, out- 
flow velocity, Doppler factor and the presence or absence 
of a re-acceleration mechanism within the hotspot. We 
consider a model in which the escape time-scale t^sc is 
the time taken for the flow to cross the hotspot and the 
cooling time-scale is the synchrotron half-life. Let R be 
the hotspot radius, v = PfC the flow velocity within the 
hotspot (note that this is not the same as the advance 
velocity of the hotspot), 5f the corresponding Doppler 
factor of the flow within the hotspot, and Ub the mag- 
netic field energy density (Ub ^ in S.I. units, 
Ub = B'^/8tt in c.g.s. units). 



2R 



'^cool ' 



7 



\dl/dt\ 
3mpC 



(11) 
(12) 

(13) 



AotU b 7 

Equating the two time-scales and combining with equa- 
tion ([8]) for the break frequency in terms of the break 
Lorentz factor, we obtain the following expression for 
the break frequency 



i/b-eO 



60 



7+4a 

(1+^) 
(1 + ^) 



niG 



-3 



iG 



-3 



R 

kpc 
R 

kpc 



GHz 



(14) 



GHz {a = (15) 



where Bsf=i is the magnetic field strength derived 
from SSC modeling under the assumption S{ — 1. 
For a Doppler factor S[ = [r(l — (3{ cos9)]^^ , the 
magnetic field strength estimated from SSC modeling 



a + 2 
» + l 



is reduced by a fa ctor of approximately 6^ """"^ (eg. 
iWorrall &: Birkinsha w 2006). Equation (fT5)) exhibits a 
strong dependence on the Doppler factor because of the 
strong dependence of the break frequency on the mag- 
netic field strength. 



Let us first consider the production of the cooling break 
in a hotspot associated with a strong relativistic nor- 
mal shock in which the post-shock velocity f3f « 0.3, 
(5t « 1, redshift z = 0.663, magnetic field B = 3 mG 
(as determined from SSC modeling in ^4.2.4p and ra- 
dius R=0.3kpc (half the geometric mean of the longest 
and shortest angular sizes of the 2.3 GHz LBA im- 
age). For such a model, the predicted break frequency 
is i^f, « 1 GHz. This is inconsistent with the lower limit 
from spectral modeling, i/f, > 500 GHz. Moreover, the 
break frequency estimated from spectral modeling is in- 
consistent with the proposed correlation between break 
frequency and equipartition magnetic field strength (eg. 
iBrunetti et "ani2003t ICheung et al.ll20b5i) which also pre- 
dicts a break frequency ~ 1 GHz. The discrepancy be- 
tween predicted break frequency and the observed lower 
limit would be alleviated if the magnetic field strength 
were < 0.15 of the value estimated from SSC modeling. 
However, if this were the case, the model SSC spectrum 
would over-predict the observed X-ray flux density. 

Let us now consider the effect of Doppler beaming on 
the observed break frequency as a possible means of re- 
solving this difficulty. Assuming a post-shock fiow veloc- 
ity /3f ~ (1 — (5f^^) and a spectral index of a ~ 0.5, a mod- 
erate Doppler factor S[ > 1.9 is sufficient to increase the 
predicted break frequency above 500 GHz, while main- 
taining agreement between the SSC model spectrum and 
the observed flux densities. 

There is also the possibility that distributed re- 
acceleration within the hots pot is affecting the produc - 
tion of the cooling break. iMeisenheimer et al.l ()1997l ) 
have suggested that distributed re-acceleration is re- 
quired to explain the spectra of the so-called low-loss 
hotspots. These are hotspots whose spectra are charac- 
terized by a power-law with a ~ 0.6 — 0.8 that extends 
to high frequency {i' > 10^^ Hz) without the predicted 
break in the spectrum. Distributed re-acceleration has 
also been proposed to explain the diffuse infra red/optical 
emission observed around some hotspots (Prieto et ciLl 
l2002tlRoeser fc Meisenheimeijll987l:|]VIeisenheimeij|2003D . 



as well as the variation in the X- ray spectral index around 
the hotspots of Cygnus A (jBalucinska-Church et al.l 
l2005f ). and the existence of flat radio spectrum re- 
gions distri buted over rnuch of the hotspot area in 
Cygnus A (jCarilli et al.l Il999( l. The favoured mech- 
anism for re-acceleration is stochastic (second or- 
der Ferm i) acceleration via magnet o hydrodynam i c tur- 
bulence (Meisenheimer et alJ Il997t iPrieto et all 120021 : 
iBal ucihska-Church et al. 200H). 

Lastly, it is possible that there is more than one site of 
particle injection, or that the hotspot is not in a steady 
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7. 



IS DOPPLER BEAMING SIGNIFICANT IN THE 
NORTHERN HOTSPOT OF PKS 1421-490? 



The aim of this section is to assess the hkehhood that 
emission from the the northern hotspot of PKS1421-490 
is Doppler beamed. To do so, in t ^7.1l we lay out the evi- 
dence for Doppler beaming in the northern hotspot, then 
discuss results of numerical simulations and radio stud- 
ies that indicate Doppler beaming may be significant in 
hotspots of radio galaxies and quasars. In H7.'2\ we con- 
sider the angle to the hue of sight of PKS1421-490, which 
should be small if Doppler beaming is to be important. 
In ^li^ we estimate the magnitude of the Doppler factor 
that would be required to account for various properties 
of the hotspot. In i j7.4l we discuss two possible arguments 
against Doppler beaming. 

7.1. Arguments for Doppler Beaming 

The northern hotspot of PKS 1421-490 is extremely 
luminous at both radio and X-ray wavelengths. The 
X-ray luminosity between 2 and 10 keV, L2-iokcV — 
3 X 10**^ ergss"^, is comparable to the X-ray luminos- 
ity of the entire jet of PKS 0637-752, without rela- 
tivistic corrections. The peak radio surface brightness 
is hundreds of times great er than that of t he bright- 
est hotspot in Cygnus A (iCarilh et al.l[l999l ). Conse- 
quently, the equipartition magnetic field strength for a 
Doppler factor of unity is greater by a factor of ~5 - 
10, and the minimum energy density is greater by a fac- 
tor of ^50 than values typically evaluated for bright 
hotspots in other radio galaxies ([M eiscnhe imer et al.l 
119971: iTavecchio et"all 120051: iKataoka fc Stawar j 12005^ 
The northern hotspot contributes 60% of the total source 
flux density at 2.3 GHz, and 75% at 8 GHz. Identifying 
the peak of region C in the ATCA image as the counter 
hotspot, we estimate the hotspot to counter hotspot flux 
density ratio to be Rhs ^ 300 at 20 GHz. In the Chandra 
X-ray band, the counter-hotspot is undetected, and we 
conservatively estimate Rhs > 100 at X-ray wavelengths. 
These are all indications that the hotspot emission may 
be Doppler beamed. Moreover, we have shown in 
that Doppler beaming may account for the absence of a 
synchrotron cooling break below 500 GHz. We now dis- 
cuss the results of numerical simulations and radio stud- 
ies that indicate Doppler beaming may be important in 
hotspots of radio galaxies. 

Numerical simulat ions of supersoni c jets in 2 and 
3 dimensions (eg. lAlov et al.l fl999l: IN orma 3 [19961 : 
iKomissarov fc Fallel 119961 : iTregillis et al.l |2001D suggest 
that flow speeds in and around hotspots can be much 
larger than those expected from the ID strong shock 
model, because the shock structure at the jet termina- 
tion is more complex than a single terminal Mach disk. 
The simulated jets undergo violent structural and veloc- 
ity changes near the jet head due to pressure variations 
in the turbulent cocoon. These violent changes in the jet 
affect the hotspot structure, and may r esult in an oblique 
shock (or sho cks) near the hotspot (jAloy et al.l 119991 : 
iNormanI 119961 ). The post-shock velocity of an oblique 
shock can be much higher than the post-shock velocity 
of a normal shock if the angle between the flow velocity 
and the shock normal is close to the Mach angle. There- 
fore, the instantaneous flow velocity through the hotspot 



may be high enough t o produce significant Doppler ef- 
fects ()Alov et al.lll999D . 

If the terminal shock is not highly oblique, the post- 
shock velocity may be relativistic if the jet contains a 
dynamically important magnetic field. The magnetic 
field can reduce the shock compression ratio and result 
in a higher post-shock Loren tz factor than that in an 
un-magnetized shock (see eg. [Double et al.l I2004D . The 
post-shock velocity in a magnetized shock depends on 
the the angle between the magnetic field and the shock 
plane, the equation of state in the pre- and post-shock 
plasma, and the mag netization parameter a = B'^/Anpc^ 
([Double et al.l 12004 ). In the case where the magnetic 
field is perpendicular to the jet direction, significant post- 
shock Lorentz factors > 2) can be achieved if ct > 3, 
depending on the equation of state. We suggest that, 
given the high magnetic field strength in the northern 
hotspot of PKS 1421-490, magnetic cushioning of the 
terminal shock due to the presence of a strong magnetic 
field in the jet may be important. 

In addition to the results of numerical simulations, ob- 
servational evidence also indicates that Doppler beam- 
ing of hotspot emission may be significant. For exam- 
ple, the brighter and more compact hotspot is gener- 
ally found on the side of the source with the brighter 
kpc-s cale jet (eg. [Bridle et al.[ [1994[ : [Hardcastle et ah] 
[19981 ). This effect is more evident in samples of quasars 
than in samples including low power sources, which 
suggests that the observed correlation between hotspot 
brightness an d jet brig htness is related to Doppler 
beaming (Har dcastle et a l. 1998). However. [Hardcastld 
([2003f ) suggest that only moderate hotspot flow veloc- 
ities {(3 ^ 0.3) ar e required to account for t his ob- 
served correlation. [Dennett-Thorpe et al.l ([19971 ) found 
that regions of high surface brightness in the lobes 
of radio galaxies have flatter radio spectra on the 
side corresponding to the brighter jet. They suggest 
that Doppler shifting of a curved hotspot spectrum 
may produce such a correlation. Again, only moder- 
ate flow speeds of 0.5 a re required to a ccount 
for this correlation ([ishw ara-Ch andra fc Saikid [2000). 
[Georganopoulos fc Kazanasi (2003) have suggested that 
deceleration of a relativistic flow from P ~ 3 to F ~ 1 in 
hotspots can explain the wide range of observed hotspot 
SEDs as bei ng pure l y an effec t of source inclination . 
However, Hardc^tld (|2003[ ) and [Hardcastle et all ([2004[ ) 
have contested this interpretation. Rather, they ar- 
gue, the shape of the hotspot SED depends only on the 
hotspot radio luminosity. 

7.2. Jet Inclination Angle 

We now consider the angle to the line of sight for 
PKS 1421-490, if Doppler beaming is to be important. 

The active galactic nucleus of PKS 1421-490 exhibits 
broad emission lines (see ij3.ip . On the basis of the 
unified scheme for active galaxies, we therefore expect 
the angle to the line of sight to be less than ^45° 
(|Urrv fc Padovarii|[19950 . Another indication of a small 
angle to the line of sight is the existence of a 60° bend 
in the northern jet, approximately 5 arcseconds (35kpc) 
from the AGN at the western end of the ridge of emis- 
sion extending west from the hotspot in Figure [T] Such a 
large jet deflection is hard to understand if it is indicative 
of the true bending angle. The well known resolution to 
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this problem is that the jet is viewed close to the line 
of sight, and the effect of projection causes a relatively 
small jet deflection to appear much larger than it actually 
is. 

7.3. Estimates of the Doppler Factor 

We now consider the magnitude of the Doppler fac- 
tor that would be required to account for the various 
observed properties. 

Let Rhs be the hotspot to counter hotspot flux density 
ratio, Pfis the bulk flow velocity in the hotspot divided by 
the speed of light, 9 the jet angle to the line of sight, and 
a the spectral index. If the two hotspots of PKS 1421- 
490 are identical, and the difference in flux density is 
purely the result of relativistic beaming, then 



Phs COS ( 



Rhs~ - 1 



(16) 



Rhs— + 1 

fee. Ide Youn'^l2002l pg. 73). 

The observed hotspot to counter-hotspot flux density 
ratio is Rhs 300 at 20.2 GHz, hence: ^hsC0s6' ~ 0.81, 
Fhs > 1.7, and 6 < 36°. A moderate Lorentz factor 
of Fhs ^ 1-7 can account for the observed hotspot flux 
density ratio. The bend in the northern jet means that 
we cannot assume the same inclination angle for the jet 
and counter-jet, so equation does not strictly apply, 
but the above calculations serve to illustrate that the 
required Lorentz factor is not large. If the jet is angled 
close to the line of sight, the real difference in inclination 
angle between jet and counter-jet may not be large. 

It should be noted that there is a difference between 
the times at which we see the two hotspots. In the case 
of PKS 1421-490 this difference is approximately (3 x 
10^/ tan6') yrs, where 6 is the angle to the line of sight. 
In equation (|16p there is an implicit assumption that the 
brightness of the hotspots remain constant over a time- 
scale of approximately 10^ — 10^ years. 

We summarize below the estimates of the Doppler 
factor required to account for various properties of the 
hotspot. 

1. A bulk Lorentz factor F > 1.7 is required to ac- 
count for the observed hotspot to counter hotspot 
flux density ratio of Rhs ~ 300. If the bend in the 
Northern jet is such that a decrease in inclination 
angle is produced, the required bulk Lorentz factor 
is lower. 

2. A Doppler factor S > 1.9 is required to account 
for the observed lower limit on the break frequency 

> 500 GHz (see 

3. A Doppler factor ^ ~ 3 is required to reduce the 
SSC model magnetic field strength to a value com- 
parable to that c alculated for other radio g alaxies 
(~ 100- 500 ^G) (iKataoka fc Stawarzl[200l . How- 
ever, such agreement is not essential since some 
variation in the radio galaxy population would be 
expected. 

7.4. Arguments Against Doppler Beaming 

In i i4.1.1l we argued that the broad emission to the west 
of the hotspot peak (region HS2) may be associated with 





Fig. 8. — 20.16GHz ATCA contours of polarized intensity with 
polarization E-vectors overlaid. The length of the polarization E- 
vectors is proportional to the fractional polarization, with 1 arc- 
second corresponding to 15% fractional polarization. The polar- 
ization vectors have been rotated by -1.8 degrees to account for 
the observed rotation measure of ~ 140 rad/m^. The main peak 
corresponds to the hotspot (regions HSl and HS2), while the sec- 
ondary peak is associated with region Jl in Figure |3] presum- 
ably associated with a jet knot. Both regions are approximately 
8% — 10% polarized. The map peak is 0.19 Jy/beam. Contour 
levels: 3.5mjy/beam X (1, 2, 4, 8, 16, 32, 64). Beam FWHM: 0.54 
X 0.36 arcsec. The scale of this image is 7.0 kpc/arcsecond. 



turbulent back-flow in the cocoon, and therefore cannot 
be Doppler beamed. As further discussed in ^ 34.1.1) if 
the high surface brightness of region HS2 is not the re- 
sult of Doppler beaming, then it seems unreasonable to 
argue that the high surface brightness of region HSl is 
the result of Doppler beaming. 

Another possible argument against Doppler beaming 
comes from interpreting the radio polarization. Figure 
[8] is a contour map of the linearly polarized intensity at 
20 GHz in the vicinity of the hotspot, with polarization 
position angle indicated by the vectors overlaid. The 
main peak in the contour map is associated with the 
hotspot (regions HSl and HS2). The offset of the sec- 
ondary peak relative to the main peak places it at the 
same position as region Jl in Figure [3] The position 
angle of the E-vectors in the hotspot indicate that the 
magnetic field (perpendicular to the E-vectors) is aligned 
nearly perpendicular to the jet direction. The magnetic 
field direction is often identified with the shock plane, be- 
cause the component of the magnetic field in the plane of 
the shock is amplified, while the component of magnetic 
field perpendicular to the shock plane is conserved. Fig- 
ure [8] therefore indicates that the terminal shock is not 
highly oblique. If the terminal shock is not highly oblique 
the post-shock velocity cannot be highly relativistic un- 
less the magnetic field is dynamically important (see §7.11 
and'Doubleiian dSOOl) ) . If the magnetic field is dynam- 
ically important, this argument against Doppler beaming 
based on the polarization position angle is not valid. 

8. INTERPRETING THE LOW FREQUENCY FLATTENING 
IN THE RADIO SPECTRUM OF THE NORTHERN 
HOTSPOT 

The aim of this section is to consider the implications 
of the observed fiattening in the hotspot radio spectrum 
discussed in §4.2.11 and illustrated in Figure [5l We con- 
sider two possible mechanisms for producing a fiatten- 
ing in the electron energy distribution at 7inin ^ 650. 
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The first mechanism we consider is the dissipation of jet 
bulk kinetic energy. Dissipation of the jet kinetic energy 
depletes low energy particles and produces a turn-over 
in the electron spectrum at a characteristic energy that 
depends on a number of parameters, including the jet 
Lorentz factor. The second mechanism we consider is a 
transition between two distinct acceleration mechanisms. 

We note that the inferred value of 7min is only weakly 
dependent on the assumed Doppler factor S because 
Doppler beaming affects the calculation of the magnetic 
field strength {B cx ^-(2+")/(i+")) as weU as the rest 
frame emission frequency corresponding to 7,nin {v' oc 
d~^i'). The value of 7min is therefore approximately pro- 
portional to 6^/^ (see Tabled]). 

8.1. Dissipation of Jet Energy 

As an illustrative calculation, we consider the dis- 
sipation of jet energy in a cold, un-magnetised pro- 
ton/electron jet. The analysis is effectively done in two 
steps: First, we use the conservation of energy and par- 
ticles to calculate the mean Lorentz factor in the hotspot 
as a function of jet Lorentz factor. We then relate the 
mean electron Lorentz factor to the peak Lorentz fac- 
tor by assuming a particular form for the electron en- 
ergy distribution. We do not specify the process by 
which the electrons and protons equilibrate. However, 
recent particle-in-cell simulations demonstrate that pro- 
tons and electrons equilibrate in un-magnetised collision- 
less shocks (Spitkovsky 2008). 

The aim of this calculation is to estimate the jet 
Lorentz factor required to produce a turn-over in the 
electron energy distribution at 7i„i,i 650 if the jet bulk 
kinetic energy is carried by protons and efficiently trans- 
fered to electrons in the hotspot. This analysis can easily 
be extended to include different jet compositions, dif- 
ferent proton to electron energy density ratios and the 
effects of the magnetic field. 

8.1.1. Model Assumptions and Definitions 

The relevant quantities are defined as follows: 7 is the 
Lorentz factor of an individual particle measured in the 
plasma rest frame, 7^ is the electron Lorentz factor at 
the peak of the electron energy distribution, F is the 
bulk Lorentz factor of the plasma, f3 — \/l ~ F"^ is the 
corresponding plasma speed in units of the speed of light 
c, 9 is the angle between the plasma velocity and the line 
of sight, S — [T{l—f3cos9)]~^ is the Doppler factor, A^(7) 
is the number density of electrons per unit Lorentz factor, 
(7) = J-fN{-f)d'y/ J N{-/)d-f is the mean Lorentz factor, 
n = J N('j)dj is the number density, e = ((7) — 1) nmc^ 
is the internal energy density, p = nm is the rest mass 
density in the plasma rest frame, p is the pressure and 
w is the relativistic enthalpy density. The relativistic 
enthalpy density is 

w = e + p + pc? (17) 

We assume the plasma is comprised of electrons (sub- 
script e) and protons (subscript p). Quantities with the 
subscript 1 refer to the jet plasma, while quantities with 
a subscript 2 refer to the hotspot plasma. 

We assume that the relativistic enthalpy density in the 
jet is dominated by the rest mass energy density of the 
the proton component, so that wi « pi^pC^. This as- 
sumption is valid provided (7)1.6 << mp/rrie. 



We assume that the electron population is ultra- 
relativistic (eg = 3pe), and that the proton population 
in the hotspot plasma is, at best, only mildly relativistic, 
and can be approximated as a thermal gas {p2,p = §e2,p)- 
We further assume that the protons and electrons equi- 
librate so that e2.p = £2,6- Hence, the hotspot pressure 
and enthalpy are p2 — £2,6 and W2 — P2(? + 3e2,e- 

8.1.2. Conservation Equations 

The equations for the conservation of energy and par- 
ticles are, respectively 

A^T\p^wi^A2TlP2W2 (18) 
AiTil3ini=A2V2l32n2 (19) 

where A is the jet cross-sectional area. Dividing equation 
(dH) by equation we can write: 
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8.1.3. Peak Lorentz Factor 



Combining equation (j20p with our model assumptions 
described above we find 



Fi«F2 1 



3e2,e 
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2 1 1 + 3^(7) 
nip 



(21) 

(22) 
and 



where we have made the substitution p2 ~ rij 

Let us introduce the parameter x= (7) / 7p which is the 
ratio of the mean electron Lorentz factor in the hotspot 
((7)) to the electron Lorentz factor at the peak of the 
electron energy distribution (7p). Then 
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In order to estimate the parameter x assume that 
the electron distribution below 7p can be approximated 
by the low energy tail of a relativistic Maxwellian, and 
above 7^ the electron distribution (before cooling via 
synchrotron emission) is a power-law extending from the 
peak to a maximum Lorentz factor 7max (see Figure |9]). 



A^(7) = 




1 < 7 < 7p 

7p < 7 < 7max 
T Tmax 



(24) 



Using this particular form for the energy distribution, 
the value of x is a function of the three parameters a, 7^ 
and 7max- If a = 2, the ratio of mean Lorentz factor to 
the peak Lorentz factor reduces to the following simple 
algebraic form 



3 3, 

h -In 

16 4 



7n 



7p 



(25) 



provided 7p >> 1. From our analysis of the spectrum of 
the northern hotspot of PKS 1421-490 in H4. 2.41 we have 
7p w 650, 7i„ax « 1-2 X 10^ and a = 2.06 so that x « 4. 
Therefore in order to produce a turn-over in the electron 
energy distribution at 7inin ~ 650 the jet must have a 
bulk Lorentz factor Fi > 5. This value of jet Lorentz 
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Fig. 9. — Plot of the assumed form of the electron energy distribu- 
tion used to relate the mean of the distribution (7) to the mode 7p 
(see equations 1124 l l and Il25l l). In i|4.2l when modeling the hotspot 
spectrum, we approximated the electron energy distribution with 
a broken power law between -f^i^ and 7max, and number density 
set to zero outside this range. Here we assume that initially the 
(injected) electron spectrum is a single (un-broken) power-law be- 
tween 7p and 7max, with N{'y) oc 7^ below 7p. The instantaneous 
cut-off at 7niin used for spectral modeling is an approximation to 
the turnover in the electron energy distribution. We identify 7niin 
with the peak of the electron energy distribution, which here we 
define as 7p. 

factor is consistent with jet Lorentz factors inferred from 
modehing t he radio to X-ray spectra of quasar jets on 
kpc-scales (iTav ecchio eFall 120001: iSchwartz et"aLl 120061: 
iKataoka fc Stawarzi,2005f ). If r2 > 1, or the jet contains 
some fraction of positron/electron pairs, or the electrons 
do not reach equihbrium with the protons, or we con- 
sider the effect of the magnetic field, then the energy 
requirements increase, and so too must Fi. 

It was noted in 21 that while only a small number of 
hotspots have provided direct estimates of 7inin, they are 
all distributed around a value of 7inin ~ 600 to within 
a factor of 2 (excluding the v alue 7 niin ~ 10* indirectly 
estimated bv lBlundell et all ()2006l )). The value of the 
parameter x is weakly dependent on the electron spec- 
trum, and in general should be within the range x 2—6. 
Therefore, dissipation of bulk kinetic energy associated 
with relativistic proton/electron jets with bulk Lorentz 
factors of order Fi > 5 can provide a natural expla- 
nation for the inferred turn-overs in electron spectra at 
7min ~ 300 - 1000. 

Our anal ysis indicates that th e value of 7min ~ 10* 
inferred bv lBlundell et all (j2006[ ) for the hotspot of the 
radio galaxy 6C0905+3955 would require a jet Lorentz 
factor Fjot > 16x + 1- If X ^ 2 (note that a ~ 0.7 for 
this hotspot, and the synchrotron spectrum extends from 
radio through to soft X-ray frequencies (jErlund et al.l 
|2008| )). then the required jet Lorentz factor is Fjet ^ 35. 

8.1.4. Electron/Positron Jet 

Let us now consider the case of a pure elec- 
tron/positron jet. In this case, assuming an ultra- 
relativistic equation of state in the jet and hotspot, the 
ratio of equations (fT5)l and implies 

Ti > (26) 

(7)l,e 

lUchivama eTal\ (|2005D estimate a mean Lorentz factor 
of (7)i,e « 50 in the jet of PKS 0637-752. If a similar 



mean Lorentz factor applies to the jet of PKS 1421-490 
then the required jet bulk Lorentz factor is Fi > 50. 

8.2. Pre- Acceleration: Cyclotron Resonant Absorption? 

We now consider an alternative explanation for the 
flattening of the electron energy distribution. The ob- 
served change in slope may be the result of a transition 
between two different acceleration mechanisms: a pre- 
acceleration process producing a relatively flat electron 
spectrum at low energy, and diffusive shock acceleration 
acting at higher energy producing an electron distribu- 
tion N{j) cx 7~" with a ~ 2. We have not modeled 
the spectrum in terms of such a scenario, but this model 
cannot yet be ruled out. 

One interesting candidate for the pre-acceler ation 
mechanism is that desc ribed bv lHoshino et al.l (|1992) and 
lAmato fc AronsI ()2006[ ). They have shown that in a rela- 
tivistic, magnetized, coUisionless shock with an electron- 
positron-proton plasma there can be efficient transfer of 
energy from protons to leptons via the resonant emis- 
sion and absorption of electromagnetic waves at high 
harmonics of the proton cyclotron frequency. This pro- 
cess produces a particle distribution described by a rel- 
ativistic Maxwell distribution at low energies (7 < Fjot) 
and a relatively flat (a < 2) power-law component ex- 
tending from ~ Fjot to ~ Fjot (wp /me). The electron 
energy index of the power-law component is sensitive 
to the plasma composition. The theoretical maximum 
Lorentz factor attained via this acceleration mechanism 
(~ Fjot('Tip/TOe)) is set by resonance with the fundamen- 
tal proton cyclotron frequency. However, the upper cut- 
off energy determined from the results of particle-in-cell 
simula tions is somewhat lowe r than the theoretical max- 
imum (jAmato fc AronsI l2006f ) . Therefore, the observed 
flattening in the hotspot radio spectrum may be asso- 
ciated with a transition between the cyclotron resonant 
absorption mechanis m and diffusive shock acceleration. 
IStawarz et al.l (|2007f ) have suggested this interpretation 
for the hotspot spectra in Cygnus A. 

9. CONCLUSIONS 

Long Baseline Array (LBA) imaging of the z=0.663 
broad line radio galaxy PKS 1421-490 has revealed 
a compact (400 pc diameter), high surface brightness 
hotspot at a projected distance of approximately 40 kpc 
from the active galactic nucleus. The isotropic X-ray lu- 
minosity of the hotspot, £2-10 kcV = 3 x 10** ergs s~^, is 
comparable to the isotropic X-ray luminosity of the en- 
tire X-ray jet of PKS 0637-752, and the peak radio sur- 
face brightness is hundreds of times greater than that of 
the brightest hotspot in Cygnus A. We successfully mod- 
eled the radio to X-ray spectral energy distribution us- 
ing a standard one-zone synchrotron self Compton model 
with a near equipartition magnetic field strength of 3 
mG. There is a strong brightness asymmetry between 
the approaching and receding hotspots, and the hot spot 
spectrum remains flat (a ~ 0.5) well beyond the pre- 
dicted cooling break for a 3 mG magnetic field, indicat- 
ing that the hotspot emission may be Doppler beamed. 
We suggest that a high plasma velocity beyond the ter- 
minal jet shock could be the result of a dynamically im- 
portant magnetic field in the jet, resulting in Doppler 
boosted hotspot emission. However, some aspects of the 
hotspot morphology may argue against an interpretation 
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involving significant Doppler beaming. LBA observa- 
tions at 1.4 GHz will be required to further investigate 
the hotspot morphology. 

There is a change in the slope of the hotspot radio 
spectrum at GHz frequencies. We successfully modeled 
this feature by incorporating a cut-off in the electron en- 
ergy distribution at 7,„in ~ 650 (assuming a Doppler 
factor of unity). If the hotspot emission is Doppler 
beamed with Doppler factor S, the low energy cut-off 
is 7min ~ 650 (5^/^. We have made use of the equa- 
tions for the conservation of energy and particles in an 
un-magnetised proton/electron jet to obtain a general 
expression that relates the peak in the hotspot electron 
energy distribution to the jet bulk Lorentz factor. We 
have shown that a sharp decrease in electron number 
density below a Lorentz factor of about 650 would arise 
from the dissipation of bulk kinetic energy in an elec- 
tron/proton jet with bulk Lorentz factor Fjet ^ 5. This 
value of jet Lorentz factor is consistent with jet Lorentz 
factors inferred from modelling t he radio to X-ray spec- 
tra of quasar j ets on kpc-scales (iTavecchio et al.ll2000l : 
ISchwartz et"all 120061 : iKataoka fc Stawarj|2005D . These 
results are of particular interest given that similar values 
of 7min have been estimated for several other hotspots. 
Our anal ysis indicates that the value of 7min ^ 10^ in- 
ferred by iBlundell et"all ()2006f ) for the hotspot of the 
radio galaxy 6C0905-I-3955 would require a jet Lorentz 
factor Fjot > 35. 

An alternative explanation for the low frequency flat- 
tening in the radio spectrum of the northern hotspot of 
PKS 1421-490 may be that it is associated with the tran- 
sition between a pre-acceleration mechanism, such as the 



cyclot ron r esonant process d escrib ed by iHoshino et al.l 
(|1992f) and lAmato fc AronsI (|2006D . and diffusive shock 
acceleration. 

Future LBA observations at 1.4GHz will help to con- 
strain the low energy end of the electron energy distri- 
bution, and infra-red observations are required to con- 
strain the high frequency end of the synchrotron spec- 
trum. More sophisticated models of the electron energy 
distribution will be required in future studies, to test the 
hypothesis that the flattening in the radio spectrum is 
associated with a transition between two distinct accel- 
eration mechanisms. 
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APPENDIX 

EQUATIONS FOR SYNCHROTRON AND SSC MODEL FLUX DENSITY 

Angle Averaged Synchrotron Flux Density for an Arbitrary Electron Energy Distribution 

The standard expressions for the synchrotron spectrum produced by an arbitrary ele ctron energy distribution con tain 
a dependence on the angle between the line of sight and the magnetic field (see eg. IWorrall fc Birkinshawll2006l ). If 
the magnetic field direction changes significantly throughout the volume in which the observed flux is produced, it is 
appropriate to use the angle averaged emission spectrum to model the source. In this appendix we present a formal 
way of calculating the angle averaged synchrotron spectrum for an arbitrary electron energy distribution. We then 
give the expression used in calculating the flux density from the volume averaged shock distribution described in the 
text. 

Let N('y) be the number density per unit Lorentz factor of relativistic electrons defined as non-zero between some 
minimum and maximum Lorentz factor 7min and 7niax in a magnetic field of flux density B. Further, define the non- 
relativistic gyro-frequency as flo, and the classical electron radius as r^- Written in S.I. units, Qq = and in c.g.s. 

units, Qq = qeB/nieC where qe and me are the electron charge and mass respectively. Let y — (47r/3)r2Q ^7"^. The 
angle averaged synchrotron emissivity (valid in both S.I. and c.g.s. units) is 



47r Ji^ 



where 



yi 



2\ 1/2 
<2 



y-^'^N{^{y))F{y)dy 

K5/3{t)dt 



(Al) 

(A2) 
(A3) 

(A4) 
(A5) 
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and the synchrotron function F{x) — x K^/^{z)dz, where Kcj/^{z) is the modified Bessel function of order 5/3. The 
integration hmits on y are given by 



47r V _2 



y2- 



7'' 

/ mm 



(A6) 
(A7) 



where and V2 are characteristic frequencies corresponding to 7niin and 7maxi viz. ^f^o7min max- 

Using the volume averaged shock distribution function N(j(y)) defined by equations ^ and Q the angle averaged 
emission spectrum is 



Dl 



A{a)iy 



-a/2 1/2 



A. 



syn / y 



'F{v)g \ y 



dy 



where 



9 y 



Aia) 




y> — 
y<- 



reTUeC 



(A8) 

(A9) 

(AlO) 

(All) 
(A12) 



(a - 1) 2°+2 

The shape of the model synchrotron spectrum is d etermined by the parameters i^i, lyb, V2 and a. The amplitude is 
governed by the parameter Asyn- In equation (jASP we have assumed the emission is produced by plasma fiowing at 

a relativistic speed /3c at an angle 9 to the line of sight with corresponding Lorentz factor V ~ {l — 0^^ and 

Doppler factor 5 = \T{\ — /3 cos throu gh a stationary volume or pattern, so that F^, cx 5^°-+'^)/'^ as appropriate for 
extragalactic jets (jLind fc Blandfordlll985l ). If the volume in which the flux is produced is moving relativistically, an 
extra factor of 6 enters, so that the leading factor of (5'"+'^)/2 in equation (|A8|1 becomes (5(°+^'/2. 

To calculate the model synchrotron spectrum, we first specify a source volume V, Doppler factor 5, redshift z and 
corresponding luminosity distance D^. We then specify lab frame values for the critical frequencies z/i, Vb and V2 
and the observation frequency v. Finally, we specify the synchrotron amplitude Asyn, and calculate the spectrum 
numerically using equation (lASp . 



Synchrotron Self Compton Flux Density 

Let Eg be the inverse Compton scattered photon energy, ei the soft photon energy, n{ei) the number density of 
soft photons per unit energy, A'^(7) the number density of relativistic electrons per unit Lorentz factor and ax the 
Thompson cross section. The inverse Compton emissivity (in the Thompson limit) from an isotropic distribution of 
relativistic electrons in an isotropic soft photon field is given by 



Scar 
167r ./o 



n (ej) 



7 



where q - 



(A13) 
(A14) 



4e,72 

and Fc(g)=-2gln(q) + q+l-2q2 

(eg. i Blumenthal & Couldl[l97l . This expression is valid for any particle distribution and any photon distribution, 
provided they are both isotropic. Changing the integration variable from 7 to q, we obtain 



Scar 
1^' 



3/2 



(ei) 



where 



e. 



91,2 



3/2 



91 



92 7 



de. 



4e,7^ ■ 

' 'mm, max 



(A15) 



(A16) 



Using the broken power law distribution described by equations ^ and ^ and expressing it in terms of the variable 
q, we find 
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jpssc 



9 • 2("-3)i? o-T 



a- 1 



-a/2 A 



(a-2)/2 



^syn 



?2 



where 



96 = 



5 - = 



1/2 



q> qt 
q<qb 



and F^f" is calculated using equation (|A8|1 . In equation (jA17|l we have assumed 

37r 



u2 rRjTv' 

h'^Vic5{l + z) 



(A17) 



(A18) 



(A19) 



(A20) 



(A21) 



as appropriate for a spherical region of homogenous plasma. In the above expression e[ and v'^ are the incident photon 
energy and frequency in the rest frame of the plasma, and j^'^ is the synchrotron emissivity at frequency v'^ in the 
rest frame of the plasma. 

We calculate the synchrotron self Compton spectrum by specifying the value of ^ssc along with the best fit values 
for Asyn, vi, Vh and 1^2 determined from fitting the synchrotron spectrum. The spectrum is calculated using equations 
(jAiyp . Together, the five parameters 1^1, Vb, 1^2, Agyn and Agsc allow the parameters K^, B, 7min, 76 and 7max to be 
determined. 
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